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I.    INTRODUCTION 

A.  BACKGROUND 

I.  Tolstoy  [Ref.  1,  Ref.  21,  and  A.  Tolstoy  and  I.  Tolstoy  [Ref.  3]  have 
recently  analyzed  the  multiple-scattering  of  sound  from  a  collection  of 
compact  (ka  «  1)  monopole  scatterers  (e.g.  gas-filled  bubbles,  gas- 
filled  balloons,  and  thin  shells)  in  water,  insonified  by  a  simple 
harmonic  plane-wave  at  frequencies  close  to  the  resonance  frequency 
of  an  individual  scatterer.  The  response  of  these  systems  was 
quantified  by  computing  the  ratio  of  the  scattered  pressure  on  the 
surface  of  a  scatterer  in  the  presence  of  the  other  scatterers  to  that  of 
the  incident  plane  wave.  This  ratio  was  defined  as  the  amplification 
factor,  \i.  They  investigated  the  value  of  this  amplification  factor  for 
various  simple  geometrical  arrangements  of  monopole  resonant 
scatterers,  and  showed  that  amplification  factors  of  the  order  of  102 
were  to  be  expected  for  air  bubbles  in  water  at  one  atmosphere 
ambient  pressure  (The  amplification  factor  of  a  single  bubble  is  about 
70.).    These  systems  were  said  to  be  "quasiresonant". 

I.  Tolstoy  [Ref.  4]  has  also  recently  analyzed  the  strength  and 
directionality  of  the  resonant  multiple-scattering  of  sound  from  a 
linear  array  of  compact  monopole  resonant  scatterers  insonified  by  a 
simple  harmonic  plane  wave.  He  reports  that  "quasiresonance"  effects 
exhibit  strong  directionality  in  this  system,  with  amplification  factors 
of  up  to  700.  The  potential  application  of  resonant  multiple- 
scattering,  or  "quasiresonance",  in  Tolstoy's  words,  to  improve  the 
performance  of  a  low  frequency  active  sonar  transducer  is  the  subject 


of  the  research  described  in  this  report. 

B.     OBJECTIVES 

The  objectives  of  this  research  are  multiple.  The  first  is  to 
reproduce  I.  Tolstoy's  results  using  network  analysis,  a  tool  familiar  to 
the  transducer  and  array  designer.  It  will  be  shown  that  network 
analysis  gives  results  identical  to  Tolstoy's  analysis.  The  existence  of 
quasiresonance  will  then  applied  to  a  circular  array  of  scatterers, 
where  the  incident  plane  wave  is  replaced  by  an  active  transducer 
located  on  the  axis  of  the  circle  and  radiating  spherical  waves.  The 
acoustic  advantage  gained  in  this  system  will  be  quantified  by  the  gain 
in  the  radiation  resistance  seen  by  the  transducer  by  the  introduction 
of  the  array  of  scatterers.  An  increase  in  radiation  resistance  enables 
the  transfer  of  more  power  to  the  acoustic  field  for  a  given  surface 
velocity  of  the  transducer  (High  power  low  frequency  active  sonar 
transducers  tend  to  be  resonant  devices  that  are  displacement 
limited.).  The  cases  of  the  transducer  in  the  plane  of  the  circle  and 
displaced  by  one-quarter  wavelength  out  of  the  plane  will  be  analyzed 
for  both  their  increase  in  radiation  resistance  and  gain  ir  directivity 
over  that  of  a  single  transducer. 
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II.    THEORY 

A.  INTRODUCTION 

For  simple  configurations  of  compact  resonant  scatterers  insonified 
by  a  simple  harmonic  plane  wave,  Tolstoy  predicted  the  existence  of 
so-called  "quasiresonance"  using  multiple-scattering  theory 
(Quasiresonance  is  Tolstoy's  term  for  the  dramatic  enhancement  of 
the  scattered  field  of  a  resonant  scatterer  as  a  result  of  resonant 
multiple-scattering.).  Tolstoy's  analysis  will  be  reviewed  and  applied 
to  the  simplest  case  of  the  scattering  of  a  plane-wave  axially  incident 
upon  a  pair  of  resonant  scatterers  in  Section  B.  In  Section  C  the  same 
problem  is  treated  using  network  analysis.  It  will  be  shown  that  this 
method  produces  results  identical  to  Tolstoy's. 

B.  TOLSTOYS  ANALYSIS  OF  RESONANT  MULTIPLE- SCATTERING 

Tolstoy's  analysis  of  the  multiple-scattering  of  a  plane  wave  incident 
upon  a  collection  of  compact  resonant  scatterers  proceeds  using  the 
linear  theory  of  acoustics.  The  following  development  parallels  that  of 
Reference  [1]  and  Reference  [3],  with  the  difference  being  that  here 
all  time-harmonic  functions  are  taken  to  be  of  the  form  d00*.  This  will 
result  in  minor  differences  in  the  appearance  of  otherwise  equivalent 
equations. 

Consider  first  the  case  of  a  single  resonant  scatterer  insonified  by 
an  incident  plane- wave.  The  scatterer  is  taken  to  be  an  air-filled 
bubble  of  radius  a,  which  is  assumed  small  compared  to  the  acoustic 
wavelength  (ka  «  1).    From  Kinsler,  Frey,  Coppens,  and  Sanders  [Ref. 
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5,  p.  230],  the  ratio  of  the  scattered  pressure  on  the  surface  of  the 
bubble  to  the  free-field  pressure  of  an  incident  plane-wave  is,  in  the 
notation  adopted  here  (note:  the  minus  sign  in  Eq.  2-1,  which  is  not 
present  in  Eq.  10-22  of  Ref.  5,  results  from  the  assumed  convention  of 
positive  velocity  directed  into  the  water): 

Plfngle(a)  = Rr  +  j(omr  (2    x) 

Pf  (Rm  +  Rr)  +  j(comr  -  s/co)  ' 

where: 

p|fngie(a)  =  scattered  acoustic  pressure  on  the  bubble  surface, 

pf  =  incident  free  field  acoustic  pressure, 

Rr  =  bubble  mechanical  radiation  resistance, 

mr  =  bubble  mechanical  radiation  mass, 

Rm  =  bubble  internal  mechanical  resistance, 

s  =  bubble  internal  mechanical  stiffness. 

The  scattered  pressure  at  any  distance  r,  plfngle^1*)'  is  related  to  the 
scattered  pressure  at  the  surface  by 

PlfngleU)  =  a  e-jkr      e-jkr       jka 
Pifngle(a)  "  r  e-Jkr  "    jkr    (1-jka)' 

where  the  approximation  ka  «  1  has  been  used.  Combining  Eq.  2-1 
and  Eq.  2-2  results  in  the  following  expression  for  the  ratio  of  the 
scattered  acoustic  pressure  at  some  radius  r  from  the  scatterer, 
p!ingie(r),  to  that  of  the  incident  pressure,  pf: 

pf      Jkrli-jd 

where: 

c  ,.(co8/cp2-i) 
S  ka 
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and 

(Oo_      /     3YPq 


CO 


pc 


Zfka)2 


(2-4) 


In  deriving  Eq.  2-3,  the  following  assumptions  and  substitutions  have 

been  made: 

(1)  Surface  tension  is  neglected,  so  the  only  mechanical 
stiffness  of  the  bubble  is  due  to  internal  gas  (air)  pressure; 

(2)  The  bubble  internal  mechanical  losses  due  to  viscosity  and 
thermal  conductivity  are  neglected  compared  to  the  acoustic 
radiation  loss,  i.e.  Rm  «  Rr; 

(3)  Rr  =  co0mr/Q; 

(4)  Q  =  l/(ka); 

(5)  s  =  (Oo2mr; 

(6)  ka  «  1; 

(7)  co  =  co0. 

It  follows  from  Eq.  2-3  that  the  ratio  of  the  scattered  acoustic 
pressure  on  the  single  scatterer's  surface  to  that  of  the  incident  free- 
field  pressure  can  be  written: 

pffngle  =  j_  i_1\  e_jka  (25) 

Pf       JkaU-Kl 

Tolstoy  introduces  the  quantities  Fa  and  A  by  rewriting  Equation  2- 
3  as  Equation  2-6. 

Fa  =  Plfngle(r)=Aedkr> 
pf  jkr 

where: 

A  =  -l— . 

i-JC 

Fa  is  the  normalized  scattered  acoustic  pressure  for  a  single  scatterer 
insonified  by  an  incident  plane-wave  and  A  is  the  single-scattering 
amplitude  coefficient. 
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Now  consider  a  system  of  N  identical  interacting  scatterers 
insonified  by  a  plane-wave.  Let  pfHr)  represent  the  scattered  acoustic 
pressure.  An  equation  for  p?f(r)/pf,  equivalent  to  that  of  Equation  2-6 
for  a  single  scatterer,  can  be  written  as  Equation  2-7  for  the  multiple- 
scattering  case: 

Fbn  =  Pm=BnS^L.  (2-7) 

Fb  is  the  normalized  scattered  acoustic  pressure  from  scatterer  n, 
insonified  by  both  the  incident  plane -wave  and  the  fields  of  the  other 
scatterers,  pd  is  the  free-field  pressure  of  the  incident  plane-wave  at 
scatterer  n,  and  Bn  is  the  multiple  scattering  amplitude  coefficient  for 
scatterer  n. 

The  response  of  a  scatterer  to  multiple  incident  pressure  fields  can 
be  evaluated  by  summing  up  the  response  to  each  individual  incident 
pressure  field  [Ref.  3].  For  a  system  of  N  scatterers,  this  results  in 
Equation  2-8  for  the  scattered  field  at  the  surface  of  the  nth  resonator: 

Bn^n=An(l+     X    Bni^^dMW^fd^n,  (2-8) 

JkSn  I  m*n  jklnm  /   Jkan 

where:  ynm  is  tne  phase  difference  in  the  incident  plane  wave 
between  resonator  m  and  n,  and  m  =  1...N.  The  subscript  notation  for 
A  has  been  added  to  allow  for  a  system  of  different-sized  scatterers. 

Tolstoy  defines  the  amplification  factor  (i  (different  from  the 
scattering  amplitude  coefficient  mentioned  previously)  to  be  the 
magnitude  of  the  ratio  of  the  scattered  acoustic  pressure  on  the 
surface  of  a  scatterer,  pftc,  to  the  free -field  acoustic  pressure  of  the 
incident  plane-wave,  p£.  With  this  definition,  Equation  2-9  follows 
from  Equation  2-7: 
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Hn  = 


PSF 


_|Bnl 


ka 


(2-9) 


n 


Tolstoy  denotes  the  amplification  factor  for  a  single  scatterer  as  u<). 
Equation  2-10,  then,  follows  from  Equation  2-6. 


W)n  = 


Pimglen 


Pn 


kan 


(2-10) 


From  Equations  2-9  and  2-10,  the  gain  in  the  amplification  factor 
due  to  multiple-scattering,  \i/\io,  is  expressed  in  terms  of  Bn/An  by 
Equation  2-11: 


P?f 


Psinglen 


=  |Bn/Anl 


(2-11) 


B.l.      APPLICATION   OF  TOLSTOYS   METHOD   TO  A   SYSTEM   OF 
TWO  RESONANT-SCATTERERS 

As  an  example  application,  Tolstoy  considered  a  system  of  two 
identical  resonators  insonified  by  a  plane-wave  traveling  along  the  axis 
joining  their  centers,  as  shown  in  Figure  2-1. 


Incoming  Plane  Wave 


Resonator  of 
radius  =  a 


© 


Figure  2-1.    A  two-resonator  system  insonified  by  a  plane-wave 
along  the  axis  of  the  resonators. 


Application  of  Equation  2-8   to  the   situation  shown   in   Figure   2-1 
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results  in  the  set  of  Equations  2-12: 


B,  £^  =  A£^  +  A^^B2  £^-  e-Jkl  , 
jka  jka  jka  jkl 


(2-12) 


B2  ed^  =  A^^  +  A^J^Bi  S=JH  e-Jkl 


jka 


jka 


jka  jkl 


For  a  system  of  Identical  scatterers,  the  values  of  Ai  and  A2  are  equal; 
therefore,  the  subscript  on  A  in  Equation  2-12  has  been  omitted. 
Since  all  variables  in  Equation  2-12  are  known  except  Bi  and  B2,  one 
can  solve  for  Bi/A  in  terms  of  A  and  kl,  resulting  in  Equation  2-13: 


Bi/A  = 


1  +(Ae-JkIje-Jkl 
jkl  1 

e-jkii2 
jkl 


(2-13) 


1-    A 


Substituting  the  expressions  for  £  and  co/coo  for  an  air-filled  bubble  into 
the  expression  for  A,  Bi/A  is  given  in  terms  of  the  variables  ka,  kl,  c, 
Po,  and  p  by  Equation  2-14: 


1  -  e-Jkl 


Psingle 


/             ka 
kl 

ej(ka- 

-kl) 

pc  ((ka)2  + 

j  ka- 

3yP0 
pc2(ka) 

)) 

\ 


/             ka 
kl 

ej(ka- 

-kl) 

pc  ((ka)2  + 

jfka- 

3yPo  \\ 

pc2(ka)  // 

(2-14) 


For  bubbles  just  under  the  surface  of  the  ocean: 

Po    =  10132  Pa  (1  atmosphere) 

p    =  998  kg  m-3 
c    =  1498  m  s"1 

As  observed  by  Tolstoy  [Ref.   3],   plots  of  Equation  2-14  display 
spurious  peaks.    Following  Tolstoy,  then,  the  data  used  to  produce  the 


16 


graphs  of  IBi/AI  versus  ka  and  kl,  plotted  in  Figures  2-2a,b,  have  been 
treated  to  limit  peaks  in  bounded  intervals  of  width  less  than  2ka 
(According  to  Tolstoy,  peaks  narrower  than  2ka,  where  2a  equals  the 
resonator  diameter,  are  not  meaningful,  since  the  pressure  over  a 
compact  resonator  is  approximately  constant.).  The  process  used 
involved  comparing  three  consecutive  values  of  IBi/AI  in  the  kl 
direction  for  each  value  of  ka  (Three  consecutive  points  in  the  kl 
direction  span  approximately  2ka.).  If  the  value  of  IBi/AI  for  the 
middle  point  was  greater  than  either  of  its  neighboring  points  (a  peak) 
then  this  value  was  set  equal  to  the  mean  of  the  two  neighboring 
points.  This  procedure  was  repeated  several  times  to  obtain  the 
desired  smoothing. 

Figures  2-2a,b  show  a  maximum  response  of  IBi/AI  equal  to  5.5  at 
approximately  ka  =  0.01393  and  kl  =  0.60.  As  a  result  of  multiple- 
scattering,  then,  the  scattered  pressure  amplitude  of  the  first  bubble 
is  5.5  times  greater  in  the  presence  of  the  second  bubble  than  in  its 
absence.  Additionally,  the  value  of  ka  at  the  maximum  response  is 
shifted  from  that  for  a  single  resonant  bubble  of  0.01379. 
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a. 


Figure  2-2a,b.  Surface  (a),  and  contour  (b)  plots  showing  IBi/AI 
for  the  double-resonator  system  as  a  function  of  ka  and  kl  for 
axial  incidence  of  the  plane-wave. 
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C.     NETWORK  ANALYSIS  OF  RESONANT  MULTIPLE-SCATTERING 

Figure  2-3  displays  a  network  model  of  N  interacting  scatterers 
insonified  by  an  incident  acoustic  wave  (Note:  pressure  (p)  and 
volume  velocity  (Un)  are  chosen  as  the  mechanical  variables,  and  the 
direction  of  positive  velocity  is  clockwise.). 


Resonators  Acoustic  Field 


Ui^ 

A 

^11 

m=  1...N 

Pi 

/ 

1 

Pi  +  X,  z]mum     \ 

u2 

V      »"i           J 

A 

^12 

P2 

r 

1 

P-,  +  L  Zom  um    A 

r\          m»2                  J 

1 
1 

i 
i 
i 

*k 

-► 

rydS 

PN 

( 

1 

PN+    ^    4LUm       ) 

K 

Figure    2-3.       Network    model    of   N    interacting    resonators 
insonified  by  a  plane -wave. 


The  network  equations  for  Figure  2-3  are: 


N 


Pn  =  Pn  +     X     ZnmUm        pn  =  -  Z&  Un  ,  (2-15) 

m=  1 


where: 

Pn    =  Surface  pressure  at  element  n, 
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pk    =  Free  field  pressure  at  element  n  due  to  an 

externally  incident  field  only, 
Un    =  Volume  velocity  of  element  n, 
Zan   =  Acoustic  impedance  of  element  n, 

Z^n   =  Acoustic  radiation  impedance  seen  by  element 

n  with  all  other  elements  blocked  (Um  *  n  =  0;    m=l...N), 

Znm   =  Acoustic  transfer  impedance  from  element  m  to 
element  n. 

Neglecting  the  effects  of  surface  tension,  thermal  relaxation,  etc., 
the  intrinsic  acoustic  impedance  of  a  small  (ka  «  1)  gas-filled  bubble 
can  be  approximated  by  one  of  reactance  only,  given  by  Kinsler,  Frey, 
Coppens  and  Sanders  [Ref.  5]  as  Equation  2-16: 

28 -'aS-^1-  (2-16' 

where: 

an  =  radius  of  element  n; 

c  =  speed  of  sound  in  water; 

k  =  wave  number,  co/c; 

P0  =  ambient  pressure  of  environment; 

Sn  =  surface  area  of  element  n; 

y  =  ratio  of  specific  heats  (air). 

The  self  acoustic  radiation  impedance  seen  by  element  n,  Z^n,  can 
be  approximated  by  Equation  2-17. 

7^n  =  pc  ((kanP+jkan)  S"1  ,  (2-17) 

where  p  is  the  density  of  water.  Equation  2-17  is  the  acoustic 
radiation  impedance  presented  to  a  compact  source  with  no  other 
active  or  scattering  elements  present.  It  is  also  the  acoustic  radiation 
impedance  seen  by  a  compact  source  in  the  presence  of  other 
compact  sources  when  these  sources  are  blocked  (Um  „  n  =  0;    m  = 
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1...N). 

The  acoustic  transfer  impedance  from  element  m  to  element  n  is 
given  by  Equation  2-18: 


Z£rm  =  Z^rm^eJ(kam  -  Unm) 
klnm 


(2-18) 


In   matrix  form,   the   network  equations   (Equation  2-15)   can  be 
written  as  Equation  2-19: 


!Pn!  =  [Z&J  (Um)  +  jpnj  ;    m,n  =  1...N 


(2-19) 


where: 


\Xn\;  n    =  1...N        A  lxN  column  matrix,  N  elements, 

[Xnm];       n,m   =  1...N        A  NxN  square  matrix,  N2  elements. 

Substituting  pn  =  -  Z^Un  and  rearranging  results  in  Equation  2-20: 


1-  Pn)  =  [Znrm  +  ^6^]  (Ura}  ;    m.n  =  1...N  , 


(2-20) 


where  8nm  is  the  Kronecker  delta  function: 


>nm  - 


1;    m  =  n  I 

\     0;    m  *n   j 


C.l.      APPLICATION   OF  NETWORK  ANALYSIS  TO  A  SYSTEM  OF 
TWO  RESONANT-SCATTERERS 

Modeling  the  system  of  two  scatterers  shown  in  Figure  2-1   as  a 
network,  Equation  2-20  becomes  Equation  2-21: 


P^ 


zH  +  z?     zf  5 
zM     zj&  +  zi 


Ui 

u2 


(2-21) 


Using  the  relation  pnc  =  ZnnUn  to  eliminate  Un,  the  above  equations 
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can  be  written  in  terms  of  the  scattered  acoustic  pressure  at  each 
individual  resonator,  pftc,  and  the  incident  plane-wave  pressures,  as 
given  by  Equation  set  2-22: 

pf  =  (p^  +  Ci2P^c)D-i1  ,  (2-22) 

P3C  =  (p£  +  C21  Pf  )  D21  , 


where 


'nm  — 


7ar 
-7a  r 


and 


rar    ,   rja 
Dn  =  - 


•yd.1       .     7a 


For  a  system  of  identical  resonators,  Cnm  and  Dn  are  the  same  for  all  m 

and   n;   therefore   the   subscript  notation   is   not  required.      Solving 

Equation  2-22  for  the  scattered  acoustic  pressure  at  each  resonator  in 

terms  of  the  incident  plane-wave  pressure  results  in  Equation  2-23  for 

the  scattered  pressure  at  resonator  one: 

,  p[  D  j-  rf  C 
D2-C2 

The  scattered  acoustic  pressure  at  the  first  resonator  in  the 
absence  of  the  second  resonator,  pffngie,,  is  obtained  using  Equation 
2-23  by  letting  the  transfer  coefficient  C  equal  zero,  resulting  in 
Equation  2-24: 

Pl&igte!  =  4  *  (2"24) 

Dividing  equation  Equation  2-23  by  Equation  2-24  results  in 
Equation  2-25  for  the  ratio  Bx/A  as  defined  by  Tolstoy: 
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Bi/A  = 


_      Pf      _ 


»sc 


Pov- 
single  i 


i+(rf/rf)(g 

i-(c" 


(2-25) 


Inserting  Equations  2-16  thru  2-18  into  Equation  2-25,  one  obtains 
an  equation  for  Bi/A  in  terms  of  the  variables  ka,  kl,  c,  Po,  pf ,  p£,  and 
P,  as  given  by  Equation  2-26: 


l-(rf/rf) 


Bi/A  = 


_      Pf      _ 


kaej(ka  -  kl) 
kl 


pc  (ka)2  +  j  ka- 


3yPo 
pc2(ka) 


ksc 


Psingle] 


kaej(ka  -  kl) 
kl 


2 


(2-26) 


pc   (ka)2  +  j  ka- 


3yP0 
pc2(ka) 


The  incident  plane  wave  pressures  p{  and  p£  are  of  equal  magnitude 
and  differ  only  by  the  phase  factor  eJ^^,  resulting  in  Equation  2-27: 


P£=    e-jkl 

pi 


(2-27) 


Substituting  this  expression  into  Equation  2-26,  one  obtains  Equation 
2-28: 


1  -  e-JW 


bi/a  =  -,££_ 

Psingle  i 


l             ka 
kl 

eJ(ka- 

-kl) 

pc  ((ka)2  + 

j(ka- 

3yPo  \\ 

pc2(ka)  // 

1  - 


ka  ^(ka  -  kl) 
kl 


\2 


(2-28) 


pc   (ka)2  +  j  ka  - 


3yPp 
pc2(ka) 


Equation    2-28    is    identical    to    Equation    2-14,    demonstrating   the 
equivalence  of  Tolstoy's  method  and  network  analysis. 
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D.     NETWORK  ANALYSIS  OF  A  TRANSDUCER  INTERACTING  WITH 
MULTIPLE  RESONANT-SCATTERERS 

Network  analysis  will  now  be  applied  to  the  problem  of  a  compact 
transducer  in  the  presence  of  an  array  of  interacting,  compact 
resonant  scatterers.  The  procedure  is  the  same  here  as  in  section  C, 
except  that  the  incident  plane  wave  is  replaced  by  an  outward- 
traveling  spherical  wave  whose  source  is  a  compact  transducer.  A 
network  model  of  the  transducer  and  N  resonator  system  is  shown  in 
Figure  2-4. 
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Figure    2-4.      Network   model   of  a    single   transducer   and   N 
interacting  resonators. 


The  network  equations  for  Figure  2-4  are: 
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po  =  TI  -  Zao  U0  ,  pn  =  -  Zn  Un  ;    n  =  1..N. 

Pn  =  X  Z^mUm  ;    m,n  =  0..N, 


(2-29a) 
(2-29b) 


m 


where 

T    =  transduction  coefficient; 
Zao   ■  transducer  open  circuit  (I  =  0)  acoustic  impedance, 

and  the  subscript  n  =  0  denotes  the  transducer,  n  *  0  a  resonator. 

The  transducer  open-circuit  acoustic  impedance  (Zao)  and 
transduction  coefficient  (T)  are  shown  only  for  completeness  and  will 
not  be  considered  further,  as  the  quantity  of  interest  will  be  the 
radiation  impedance  presented  to  the  transducer. 

Using  the  network  Equations  2-29,  the  matrix  equation  for  the 
network  model  shown  in  Figure  2-4  is  given  by  Equation  2-30: 


PO 

(Pn! 


m  (zgfo!T 

(Znb|     L^nmJ 


u0 

!U 


mi 


n,  m  =  1...N 


(2-30) 


where 


IXnl,  n  =  1...N  denotes  a  lxN  column  matrix,  N  elements, 
[Xnm],  n,m  =  1...N  denotes  a  NxN  square  matrix,  N2  elements, 

and  the  superscript  T  denotes  transpose.  Replacing  the  1  x  N  surface 
pressure  matrix,  {pnK  with  {-  ZnUn},  and  then  normalizing  by  Uo,  the 
transducer  volume  velocity,  Equation  2-30  becomes  Equation  2-31: 


Zaor 


tal 


10! 


where 


m 

(Zno! 


\ZmOi 
Znm  +  ZnOntnJ 


lTTreH 


;  n,m  =  l...N,(2-31) 


Zfotal    =  total  acoustic  radiation  impedance  seen  by  the 
transducer; 
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Um1   =  ratio  of  element  m's  volume  velocity  to  that  of  the 
transducer  (U^1  =  Um/U0). 

Equation  2-31  can  be  written  as  two  separate  matrix  equations, 
Equations  2-32a,b: 

Zf0rtai  =  m  +  (zSbHug1!  •  (2-32a) 

{0}  =  (Z&bj  +  [Zgm  +  ZftfnnJ  (u£!i  .  (2-32b) 

Solving  Equation  2-32b  for  |Uml  by  matrix  inversion  and  substituting 
the  result  into  Equation  2-32a,  one  obtains  Equation  2-33,  the  total 
acoustic  radiation  impedance  seen  by  the  transducer. 

ZfJtal  =  Z6S  "   (ZgS>|T  [zfSn  +  Z*8n  J"1  (Zgfc)  .  (2-33) 

Equation  2-33  expresses  this  impedance  in  terms  of  the  radiation  and 
transfer  impedances  of  the  individual  elements. 

The  total  acoustic  power  radiated  by  the  transducer  is  given  by 
Equation  2-34. 

n  =  i|Ud2Re(Zf0rtai).  (2-34) 

The  gain  in  power  delivered  to  the  acoustic  field  due  to  the  addition  of 
the  sympathetic  resonators  equals  the  ratio  of  total  power  delivered 
with  the  resonators  present  to  the  power  delivered  when  they  are 
absent,  IIo  =  —  |  Uo|2  Re(Z§5  )•  For  a  displacement-limited  source, 
typical  of  a  sonar  projector  operated  near  its  resonant  frequency,  this 
ratio  reduces  to  the  ratio  of  the  radiation  resistances  for  the  two  cases. 
This  gain  factor  is  expressed  as  Equation  2-35: 

Gr  "1^85r  ,2"35) 
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III.   ANALYSIS 

A.  INTRODUCTION 

The  theory  developed  in  the  previous  chapter  will  now  be  applied 
to  two  specific  geometric  configurations  of  a  transducer  and  N 
identical  resonators:  a  "planar"  and  a  "conical"  configuration.  The 
resonators  will  be  treated  as  simple  gas-filled  bubbles.  It  is  assumed 
throughout  that  the  acoustic  wavelength  is  much  larger  than  the 
dimensions  of  either  the  transducer  or  the  resonators.  For  initial 
calculations,  ambient  pressure  is  taken  to  be  one  atmosphere,  and  the 
radius  of  the  transducer  is  taken  to  be  equal  to  that  of  the  resonators 
(kao  =  kan).  Later,  calculations  will  be  presented  for  the  case  of  an 
ambient  pressure  of  50  atmospheres,  and  the  restriction  kao  =  kan  will 
be  relaxed. 

B.  ANALYSIS  OF  PLANAR  CONFIGURATION 

The  geometry  of  the  first  arrangement  to  be  analyzed  is  shown  in 
Figure  3- 1 .    It  consists  of  N  identical  resonators,  equally  spaced 
around  a  circle  of  radius  lno-     This  geometric  arrangement  will  be 
referred  to  as  the  planar  configuration. 

For  the  configuration  shown  in  Figure  3-1,  Equation  2-35  becomes 
Equation  3-1: 

q     _  Re\  Zqq  -    )Zmo|     [Znm  +  ZnOnniJ      |ZnO|J  •    n  m  =  1     4  (3    1) 

~Re(Zg5)  -•      ■     - 
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Figure  3-1.    Resonators  equally  spaced  around  transducer 
in  a  planar  configuration  of  radius  lno.  N=4. 


Using  Equations  2-16,  17,  and  18,  and  canceling  the  common  surface 
areas  Sn,  Equation  3-1  can  be  written  in  terms  of  specific  acoustic 
impedances  as  Equation  3-2: 

_  Re(za5  -   (zMo|T  [zftfa  +  zftSnm]"1  (zg6i)  .    n  m  .,  1      A  fo   91 

UR~~  Relzft)  '    n'm"  1""4  [S'2) 

where: 

z55  =  pc  ((kaj2  +jka), 

zftm  +  zftSnm  =  pc  ({kaj2  +  jka)  -  j    J  °  ;    n  =  m, 

c(ka) 

Zftfn  +  Zft5nm  =  Zftfn  =  pc  ((ka)2  +  jka)  -£*-   d  (^  "  klnm)  ;    n  *  m, 

klnm 

zgfc  -  pc  ((kaj2  +  jka)  -^a_   d  (ka  -  kln0)      , 

klnO 

and  lnm  is  the  distance  between  element  n  and  m,  the  subscript  0 
denoting  the  transducer.  Using  simple  trigonometry,  klnm  for  n,m  * 
0  can  be  written  in  terms  of  m,  n,  and  kln0  for  the  planar  case  as 
Equation  3-3: 


28 


klnm  =   kln0  2  sinj-E-  |n  -  m|)  ;    n,m  =  1...N.  (3-3) 

Using  Equations  3-2  and  3-3,  the  gain  in  radiation  resistance  seen  by 
the  transducer,  Gr,  can  be  expressed  in  terms  of  ka  and  klno- 

The  gain  in  radiation  resistance  for  the  N  =  4  planar  configuration 
can  be  plotted  in  the  form  of  surface  and  contour  plots  using  Equation 
3-1,  as  shown  in  Figures  3-2a  and  2b.  It  should  be  pointed  out  that 
here  and  in  all  subsequent  surface  plots  of  Gr,  no  spurious  peaks  were 
observed,  and  no  smoothing  has  been  applied.  Inspection  of  these 
graphs  reveals  two  local  maxima,  one  of  magnitude  4.34  at  ka  =  0.975 
kar,  klno  =  0.5,  and  another  of  magnitude  1.79  at  ka  =  1.005  kar,  klno 
=  2.7,  where  kar  equals  0.01379,  the  resonance  value  of  ka  for  a  single 
air-filled  bubble  at  one  atmosphere  ambient  pressure.  The  local 
maximum  at  ka  =  0.975  kar,  klno  =  0.5  is  an  artifact  of  the  range  of 
klno  chosen.  In  fact,  the  gain  in  radiation  resistance  diverges  as 
(klno)  as  klno  approaches  zero.  This  case  is  uninteresting  and  will 
not  be  considered  further.  The  interesting  local  maximum  at  ka  = 
1.005  kar,  klno  =  2.7  represents  a  gain  of  almost  a  factor  of  two  in  the 
radiation  resistance  seen  by  the  transducer.  Note  that  kl  =  n 
corresponds  to  a  separation  distance  of  one-half  wavelength.  For  the 
N  =  4  planar  configuration,  then,  the  maximum  gain  in  radiation 
resistance  is  obtained  when  the  resonators  are  located  a  distance  just 
under  one-half  wavelength  from  the  transducer. 

The  gain  in  radiation  resistance  for  the  planar  configuration  where 
N  is  set  equal  to  eight  and  fifteen  is  shown  in  Figures  3-2c  thru  3-2f. 
These  plots  display  similar  characteristics  to  the  N  =  4  case.  All  show 
a  local  maximum  of  magnitude  approximately  two  at  ka  =  kar  = 
0.01379  and  kln0  =  n. 
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The  local  maximum  values  for  Gr  at  klno  =  n,  ka  =  kar  =  0.01379  for 
various  values  of  N  ranging  from  2  to  25  and  their  associated 
parameters  are  tabulated  in  Table  3-1.  Two  major  points  can  be 
observed  from  Table  3-1.  First,  the  resonant  bubbles'  volume  velocity 
leads  the  transducer's  by  a  phase  shift  of  90°  (equivalent  to  -  270°). 
This  phase  shift  is  created  in  two  parts.  Since  the  separation  between 
transducer  and  resonator  is  approximately  one-half  wavelength,  this 
introduces  a  shift  of  -  180°  in  phase  due  to  simple  wave  propagation. 
The  remaining  -  90°  is  introduced  at  the  surface  of  a  compact 
resonating  bubble  between  the  incident  and  scattered  wave,  as 
predicted  by  Equation  2-24.  Second,  note  that  little  is  gained  as  one 
increases  the  number  of  resonators  beyond  about  four. 


N 

Gr 

K^max 

J^Mnax 

magOJJf1) 

aigOJJf1) 

2 

1.26 

0.01382 

2.65 

0.395 

90.0° 

4 

1.87 

0.01383 

2.82 

0.649 

90.0° 

6 

1.96 

0.01385 

2.87 

0.474 

90.0° 

8 

1.96 

0.01386 

2.87 

0.356 

90.0° 

10 

1.96 

0.01386 

2.87 

0.285 

90.0° 

15 

1.96 

0.01385 

2.87 

0.190 

90.0° 

20 

1.96 

0.01381 

2.87 

0.142 

90.0° 

25 

1.96 

0.01376 

2.87 

0.114 

90.0° 

Table  3-1.  Values  of  Gr  for  N  ranging  from  2  to  25  and  the 
parameters  associated  with  the  local  maximum  at  klno  =  n, 
ka  =  kar  =  0.01379.    Planar  configuration. 
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Figure  3-2a,b.  Surface  and  contour  plots  of  Gr  as  a 
function  of  ka  and  klno  for  N  =  4.  Planar  configuration.  P0 
=  1  atm. 
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c. 


Figure  3-2c,d.  Surface  and  contour  plots  of  GR  as  a 
function  of  ka  and  klno  for  N  =  8.  Planar  configuration.  P0 
=  1  atm. 
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Figure  3-2e,f.  Surface  and  contour  plots  of  GR  as  a 
function  of  ka  and  klno  for  N  =  15.  Planar  configuration. 
P0  =  1  atm. 
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B.l.     PLANAR  CONFIGURATION  BEAM  PATTERNS 

Polar  plots  of  the  beam  patterns  can  be  developed  for  the  planar 
configuration  using  the  axis  system  defined  in  Figure  3-3.  The  origin 
is  located  at  the  transducer,  with  the  resonators  located  in  the  x,y 
plane  and  with  9  being  the  angle  measured  from  the  positive  y  axis  in 
the  x,y  plane. 

Figure  3-4  clearly  demonstrates  the  acoustic  advantage  of  the 
planar  configuration  over  that  of  a  lone  transducer.  The  zero  dB 
reference  is  taken  as  the  field  of  an  identical  transducer  with  equal 
volume  velocity,  in  the  absence  of  the  resonators. 


z  axis 


x  axis 


y  axis 


y.z  plane 


^5 


z,9  =45°plane 


Figure  3-3.  Axis  system  for  beam  patterns  plots.  Planar 
configuration.  Note  the  definitions  of  the  y,z  and  z,0  =  45° 
planes. 
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Figure  3-4.    Beam  patterns  in  the  y,z  and  z,0  =  45°  planes. 
Planar  configuration.    P0  =  latm.    Scale  in  decibels  (dB). 


C.     ANALYSIS  OF  CONICAL  CONFIGURATION 

The  second  geometric  arrangement  analyzed  was  chosen  to  take 
advantage  of  the  90°  phase  shift  between  the  volume  velocity  of  the 
resonators  and  the  transducer  for  maximum  radiated  power.  If  the 
transducer  is  offset  from  the  plane  of  the  resonators  by  one-quarter 
wavelength  along  the  positive  z  axis,  the  radiation  from  the  resonators 
will  constructively  interfere  with  that  of  the  transducer  along  the 
positive  z  axis.  This  configuration  will  increase  the  directivity  of  the 
system.  This  arrangement  will  be  referred  to  as  the  conical 
configuration,  and  is  shown  in  Figure  3-5. 
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Transducer 


Resonator  n 

Figure  3-5.  Resonators  equally  spaced  around  transducer 
in  a  conical  configuration  at  radius  R.  N  =  4.  The 
transducer  is  displaced  out  of  the  plane  of  the  resonators 

by  a  distance  of  X/4. 

The  value  of  Gr  can  again  be  expressed  in  terms  of  ka  and  kln0 
using  Equation  3-1,  where  Equation  3-4  is  used  to  express  klnm  in 
terms  of  n,  m,  and  klno  for  the  conical  configuration. 


klnm  =  ytklno)2-^)2  2sin(^|n-m|);    n,m  =  1...N.  (3-4) 

The  gain  in  radiation  resistance  for  the  conical  configuration  for  the 
cases  N  =  4,  8,  and  15  are  plotted  in  the  form  of  surface  and  contour 
plots,  as  shown  in  Figures  3-6a  through  6f.  Inspection  of  these  plots 
reveals  the  existence  of  the  same  local  maximum  as  the  planar 
configuration  of  magnitude  slightly  less  than  2  at  klno  =  n  and  ka  = 
0.01379.  Note  that,  unlike  the  planar  case,  the  kl  axis  of  Figures  3-6a 
through  f  begins  at  kln0  =  1.75  vice  0.5.  This  is  due  to  the 
displacement  of  the  transducer  in  the  positive  z  direction.  This 
displacement  requires  kln0  be  at  least  n/2. 

The  local  maximum  values  for  Gr  at  klno  s  n,  ka  s  0.01379  for 
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various  values  of  N  ranging  from  2  to  25  and  their  associated 
parameters  are  tabulated  in  Table  3-2  for  the  conical  configuration.  In 
each  case  the  value  of  Gr  for  the  conical  configuration  is  slightly  less 
than  for  the  planar  configuration.  Also,  the  value  of  klmax  is  slightly 
greater  for  the  the  conical  configuration,  but  is  still  approximately 
equal  to  n  (A  value  of  n  for  klno  corresponds  to  a  value  of  30  degrees 
for  the  angle  a  in  Figure  3-5.). 


N 

Gr 

K'^max 

^max 

magCUJf1) 

argOJ-1) 

2 

1.29 

0.01382 

2.77 

0.429 

90.0° 

4 

1.78 

0.01383 

3.08 

0.601 

90.0° 

6 

1.82 

0.01385 

3.12 

0.428 

90.0° 

8 

1.82 

0.01386 

3.12 

0.321 

90.0° 

10 

1.82 

0.01385 

3.12 

0.257 

90.0° 

15 

1.82 

0.01383 

3.12 

0.171 

90.0° 

20 

1.82 

0.01378 

3.12 

0.128 

90.0° 

25 

1.82 

0.01372 

3.12 

0.103 

90.0° 

Table  3-2.  Values  of  Gr  for  N  ranging  from  2  to  25  and  the 
parameters  associated  with  the  local  maximum  at  kl  =  n, 
ka  =  klr  =  0.01379.    Conical  configuration. 
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a. 


Figure  3-6a,b.  Surface  and  contour  plots  of  Gr  as  a 
function  of  ka  and  klno  for  N  =  4.  Conical  configuration. 
P0  =  1  atm. 
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kl 


Figure  3-6c,d.  Surface  and  contour  plots  of  Gr  as  a 
function  of  ka  and  klno  for  N  =  8.  Conical  configuration. 
P0  =  1  atm. 
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Figure  3-6e,f.  Surface  and  contour  plots  of  Gr  as  a 
function  of  ka  and  kln0  for  N  =  15.  Conical  configuration. 
P0  =  1  atm. 
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C.l.     CONICAL  CONFIGURATION  BEAM  PATTERNS 

To  demonstrate  the  improved  directivity  of  the  conical 
configuration  over  that  of  the  planar  system,  polar  plots  can  be 
developed  using  the  axis  system  as  defined  in  Figure  3-7,  where  the 
transducer  is  displaced  in  the  positive  z  direction  from  the  plane  of 
the  resonators. 


z  axis 


Figure  3-7.    Axis  system  for  beam  pattern  plots.     Conical 
configuration. 


Figure  3 -8a  shows  the  beam  pattern  of  the  conical  configuration  in 
both  the  y,z  and  the  z,0  =  45°  planes.  As  expected,  the  symmetry  in 
the  y,z  plane  about  the  z  axis  is  lost,  with  the  sound  pressure  level 
(SPL)  in  the  positive  z  direction  being  8  dB  greater  than  in  the 
negative  z  direction.  If  one  compares  the  beam  patterns  of  the  conical 
configuration  to  that  of  the  planar  configuration,  shown  in  Figure  3 -8b, 
one  observes  that  the  SPL  in  the  positive  z  direction  for  the  conical 
configuration  is  1  dB  greater  than  for  the  planar  configuration,  and  10 
dB  greater  than  for  one  transducer  alone.  Figure  3 -8b  clearly  shows 
that  significant  directional  discrimination  can  be  achieved  with  only  a 
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minor  decrease  in  total  power  gain,  Gr. 


z  axis 


0  xy  plane 


LEGEND 
Z,Y  PLHNE 


Figure  3-8a.    Beam  patterns  in  the  y,z  and  z,G  =  45°  planes. 
Conical  configuration.    P0  =  1  atm.    Scale  in  decibels  (dB). 
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Figure  3-8b.  Beam-patterns  in  the  y,z  plane  for  the  planar 
and  conical  configurations.  P0  =  1  atm.  Scale  in  decibels 
(dB). 


D.     OPTIMUM  CONFIGURATION 

The  results  of  all  calculations  performed  have  been  presented  in 
terms  of  the  dimensionless  radius  and  spacing,  ka  and  klno-  The 
actual  dimensions  for  the  radius  a  and  the  separation  distance  lmax  for 
a  given  frequency  are  obtained  using  Equation  3-5: 


-  kac       1   _  _  klnpc 
"  2nf        ln°  "  ~2kT  * 


(3-5) 


Taking  ka  =  kar  =  0.01379  and  kln0  s  n  in  Equation  3-5  for  the 
optimum  values  of  a  typical  system  at  an  ambient  pressure  of  one 
atmosphere,  a  plot  of  the  resulting  dimensions  a  and  lno  as  a  function 
of  frequency  can  be  made,  as  shown  in  Figures  3-9a  and  3-9b: 
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Figure  3 -9a.    Optimum  resonator  diameter  (a)  at  resonance 
versus  frequency.    Pq  =  1  atm. 
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Figure  3-9b.  Optimum  resonator-transducer  separation 
distance  (lno)  for  maximum  power  gain  versus  frequency. 
PO  =  1  atm. 
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E.     ANALYSES  WITH  TRANSDUCER  ka  =  CONSTANT  AND  P0  =  50 
ATM 

Two  major  assumptions  that  have  been  made  throughout  the 
analysis  up  this  point  are  now  modified.  The  first  is  that  the  ambient 
pressure  is  one  atmosphere  and  the  second  is  that  the  radius  of  the 
transducer  equals  that  of  the  resonators  (kao  =  kan).  In  the  following 
analysis,  an  ambient  pressure  of  50  atmospheres  (corresponding  to  a 
depth  of  500  meters)  will  be  taken  as  a  more  realistic  value  for  an 
operational  system.  Since  the  internal  acoustic  impedance  of  a  bubble 
is  a  function  of  depth,  the  value  for  kar  of  a  single  bubble  at  resonance 
will  shift  from  0.01379  (P0  =  1  atm)  to  0.09754  (P0  =  50  atm).  Also, 
the  value  of  kao  wiU  no  longer  be  considered  equal  to  that  of  the 
resonators,  but  will  be  held  constant  at  various  specific  values. 

In  the  previous  sections  it  was  shown  that  little  is  gained  when  the 
number  of  sympathetic  resonators  is  greater  than  about  four;  hence 
the  following  analysis  will  be  restricted  to  the  case  N  =  4.  Conducting 
similar  calculations  as  before,  contour  plots  and  surface  plots  have 
been  generated  for  the  cases  N  =  4,  kao  =  0.05,  0.10  and  0.20,  for 
both  the  conical  configuration  (Figures  3- 10a  through  3-10f)  and  the 
planar  configuration  (Figures  3- 11a  through  3-1  If).  As  in  the  previous 
analyses  where  the  ambient  pressure  was  taken  to  be  one  atmosphere, 
all  cases  display  a  local  maximum  of  GR  of  slightly  less  than  2  for  ka  = 
kar  =  0.09754  and  kl  =  n.  These  maximum  values  are  tabulated  in 
Table  3-3.  Also  tabulated,  though  not  plotted,  are  the  cases  kao  = 
0.025  and  0.150.  The  first  entry  in  Table  3-3  for  either  configuration 
was  calculated  by  holding  kao  equal  to  kan,  as  in  previous  calculations, 
vice  holding  kao  constant.    A  review  of  Table  3-3  and  Figures  3-11  and 
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12  demonstrates  that  the  value  of  kao  has  little  effect  on  the  achievable 
power  gain. 


configuration 

kao 

Gr 

^<*max 

^max 

magOJn61) 

arg(Unel) 

planar 

kan 

1.87 

0.09987 

2.82 

0.649 

90.0° 

.025 

1.87 

0.09987 

2.82 

0.652 

90.0° 

.050 

1.87 

0.09987 

2.82 

0.652 

90.0° 

.100 

1.87 

0.09987 

2.82 

0.649 

90.0° 

.150 

1.87 

0.09987 

2.82 

0.645 

90.0° 

.200 

1.87 

0.09998 

2.82 

0.639 

90.0° 

conical 

ka,. 

1.78 

0.09990 

3.08 

0.601 

90.0° 

.025 

1.78 

0.09990 

3.08 

0.603 

90.0° 

.050 

1.78 

0.09990 

3.08 

0.603 

90.0° 

.100 

1.78 

0.09990 

3.08 

0.601 

90.0° 

.150 

1.78 

0.09990 

3.08 

0.597 

90.0° 

.200 

1.78 

0.09990 

3.08 

0.592 

90.0° 

Table  3-3.  Values  of  Gr  for  various  values  of  kao  ranging 
from  0.025  to  0.2  and  the  parameters  associated  with  the 
local  maximum  at  kl  =  n,  ka  =  kar  =  0.09754.  Conical  and 
planar  configurations.    Pq  =  50  atm. 
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Figure  3-10a,b.  Surface  and  contour  plots  of  Gr  as  a 
function  of  ka  and  klno  for  N  =  4.  Planar  configuration.  P0 
=  50  atm.    ka0=  0.05. 
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c. 


Figure  3-10c,d.  Surface  and  contour  plots  of  Gr  as  a 
function  of  ka  and  klno  for  N  =  4.  Planar  configuration.  P0 
=  50  atm.    kao  =  0.10. 
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c. 


Figure  3-10e,f.  Surface  and  contour  plots  of  GR  as  a 
function  of  ka  and  kln0  for  N  =  4.  Planar  configuration.  P0 
=  50  atm.    kao  =  0.20. 
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kl 


Figure  3-lla,b.  Surface  and  contour  plots  of  Gr  as  a 
function  of  ka  and  kln0  for  N  =  4.  Conical  configuration. 
P0  =  50  atm.    ka0  =  0.05. 
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d. 


ki 


Figure  3-1  led.  Surface  and  contour  plots  of  Gr  as  a 
function  of  ka  and  klno  for  N  =  4.  Conical  configuration. 
P0  =  50  atm.    ka0  =  0.10. 
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Figure  3-lle,f.  Surface  and  contour  plots  of  Gr  as  a 
function  of  ka  and  kln0  for  N  =  4.  Conical  configuration. 
P0  =  50  atm.    kao  =  0.20. 
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E.l.     BEAM  PATTERNS,  P0  =  50  ATM 

Since  there  is  no  significant  variation  in  gain  with  the  value  of  ka0, 
beam  patterns  were  only  developed  for  the  case  ka0  =  kar,  and  are 
shown  in  Figures  3-12  through  3-14.  These  plots  show  the  maximum 
SPL  in  the  positive  z  direction  for  the  conical  configuration  to  be  1.5 
dB  greater  than  for  the  planar  case  and  10.5  dB  greater  than  for  the 
transducer  in  the  absence  of  the  resonators. 


z  axis 


xy  plane 


Figure  3-12.  Beam  patterns  in  the  y,z  plane  for  the  planar 
and  conical  configurations.  Pq  =  50  atm.  Scale  in  decibels 
(dB). 
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Figure  3-13.  Beam  patterns  in  the  y,z  and  z,9  =  45° 
planes.  Planar  configuration.  P0  =  50  atm.  Scale  in 
decibels  (dB). 
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Figure  3-14.  Beam  patterns  in  the  y,z  and  z,9  =  45° 
planes.  Conical  configuration.  P0  =  50  atm.  Scale  in 
decibels  (dB). 
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E.2.     OPTIMUM  CONFIGURATION,  P0  =  50  ATM 

The  actual  dimensions  of  the  resonator  radius  a  and  the  separation 
distances  lnQ,  for  a  system  operated  at  an  ambient  pressure  of  50 
atmospheres,  were  calculated  using  Equation  3-5  using  the  optimum 
values  ka  =  kar  =  0.09754  and  klno  =  n.  Since  the  values  of  lno  are  not 
depth-dependent,  Figure  3-9b  still  applies.  Figure  3-15  displays  the 
relation  between  the  optimum  resonator  radius  and  frequency  for  an 
ambient  pressure  of  50  atmospheres. 
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Figure    3-15.        Optimum    resonator    diameter    (a)    at 
resonance  verses  frequency.    P0  =  50  atm. 
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IV.    SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  SUMMARY 

The  problem  of  multiple -scattering  of  an  incident  plane-wave  from 
a  system  of  compact  resonant-scatterers  was  solved  using  both 
Tolstoy's  method  and  network  analysis.  It  was  shown  that  these  two 
methods  of  analysis  are  equivalent  and  that  both  predict  the  existence 
of  "quasiresonance".  Network  analysis  was  then  applied  to  a  circular 
array  of  scatterers,  where  the  incident  plane-wave  was  replaced  by 
spherical  waves  radiated  from  an  active  transducer,  located  on  the 
center  axis.  The  acoustic  advantage  of  this  system  was  quantified  by 
the  gain  in  the  radiation  resistance  seen  by  the  transducer  in  the 
presence  of  these  scatterers  compared  to  a  lone  transducer.  Two 
configurations  were  analyzed:  (1)  the  transducer  in  the  plane  of  the 
resonators,  and  (2)  the  transducer  displaced  by  one-quarter 
wavelength  from  the  plane.  Both  configurations  showed  a  gain  in 
radiation  resistance  of  almost  a  factor  of  two  and  an  improved 
directivity  over  that  of  a  single  transducer. 

B.  CONCLUSIONS 

The  numerical  analysis  of  a  transducer  surrounded  by  a  circular 
array  of  sympathetic  resonators  indicates  that  a  gain  in  radiation 
resistance  of  about  two  is  obtainable  compared  to  a  lone  transducer. 
The  physical  dimensions  of  such  a  system  operating  at  an  ambient 
pressure  of  50  atmospheres  are  displayed  in  Table  4-1  for  several 
frequencies  of  interest. 

57 


Frequency 
(Hz) 

Resonator  Radius  (a) 
(m) 

System  Radius  (R 
(m) 

conical          planar 

10 

100 

1000 

0.329 

0.0329 

0.00329 

64.9 
6.49 
0.649 

74.9 
7.49 
0.749 

Table  4-1.  Dimensions  of  planar  and  conical 
configurations  for  various  frequencies.  Po  =  50  atm.  R  is 
the  radius  of  the  array. 


Although  the  calculations  presented  here  were  for  spherical 
bubbles,  the  results  also  apply  to  more  realistic  devices  (e.g.  thin 
shelled  cavities,  balloons,  etc)  as  long  as  they  are  compact  (ka  «  1). 
M.  Strasberg  [Ref.  6]  showed  that  the  shape  of  a  resonant  body  has 
little  effect  on  its  fundamental  resonance  frequency.  For  example,  he 
calculated  that  an  oblate  spheroid  with  a  ratio  of  lengths  of  major  to 
minor  axes  of  four,  the  fundamental  resonance  frequency  differs  by 
only  a  factor  of  1.08  from  that  of  a  sphere  with  the  same  volume. 
Therefore,  any  deformation  of  the  resonators  from  a  true  spherical 
shape,  for  example  due  to  forces  such  as  buoyancy,  should  have  little 
affect  on  system  performance. 

C.      RECOMMENDATIONS 

The  results  of  the  research  reported  in  this  thesis  are  very 
encouraging.  It  is  recommended  that  follow-on  work  include  the 
design,  construction,  and  testing  of  an  array  to  verify  the  theoretical 
predictions  presented. 

Theoretical  analysis  remains  to  be  conducted  in  at  least  four  areas. 
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First,  in  all  of  the  calculations  performed,  the  sympathetic  resonators 
were  assumed  to  be  gas- filled  air  bubbles  executing  radial  oscillations 
only.  A  more  realistic  model  of  a  sympathetic  resonator  may  be 
required,  based  upon  the  properties  of  a  practical  device.  It  is 
recommended  this  be  investigated.  Second,  Tolstoy  has  shown  that 
the  scattering  of  an  incident  plane-wave  from  a  linear  array  of  resonant 
scatterers  exhibits  strong  directionality.  Perhaps  a  linear  array  of 
sympathetic  resonators  offers  better  performance  gain  than  the 
circular  arrays  considered  in  this  research.  It  is  recommended  that 
such  a  configuration  be  investigated.  Third,  it  was  shown  that  in  the 
optimum  conical  configuration  the  ring  of  sympathetic  resonators  lies 
on  a  cone  with  a  apex  angle  of  approximately  120  degrees, 
independent  of  frequency.  This  suggests  that  a  broad-band,  constant 
directivity  system  could  be  constructed  using  a  series  of  concentric 
rings  of  resonators,  each  tuned  to  a  slightly  different  frequency  and 
placed  in  its  appropriate  location  along  the  surface  of  such  a  cone,  as 
indicated  in  Figure  4-1.  It  is  recommended  that  the  properties  of 
such  a  system  be  investigated.  Fourth,  in  the  calculations  presented, 
the  sympathetic  resonators  were  all  assumed  to  be  identical,  and  their 
locations  were  assumed  to  be  precisely  known,  as  well  as  that  of  the 
transducer.  It  is  recommended  that  the  effects  of  resonator 
nonuniform! ty  and  imprecise  location  be  investigated. 


59 


Transducer 


Resonator  r=a' 


Resonator  r=a 


Figure  4-1.  Multiple-ring  conical  configuration  allowing  an 
expanded  frequency  response.  Two  concentric  rings  of 
resonators  of  radii  a  and  a',  separated  from  the  transducer 
by  distances  of  1  and  1\  respectively. 
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LIST  OF  SYMBOLS 

an  Radius  of  element  n. 

c  Speed  of  sound  in  water. 

Gr  The  ratio  of  the  acoustic  radiation  resistance  of  a  transducer 

with  sympathetic  resonators  present  to  the  acoustic  radiation 
resistance  of  an  identical  transducer  with  no  resonators 
present. 

k  Wave  number,  co/c. 

kar  The  resonance  value  of  ka  for  a  single  air- filled  bubble. 

Inm  Separation  distance  between  element  n  and  m. 

pfx  Incident  free  field  acoustic  pressure  at  element  n  due  to 

external  field  only. 

pn  Surface  pressure  at  element  n. 

pftc  Scattered  pressure    at  the  surface  of  resonator  n,  in  the 

presence  of  other  resonators. 

Pimgle     Scattered  pressure  at  the  surface  of  a  single    resonator,  i.e.  in 
the  absence  of  other  resonators. 

r  Distance  from  scatterer  center  to  a  field  point. 

z  Specific  acoustic  impedance,  pressure  divided  by  surface 

normal  velocity. 

A  Single-scattering  amplification  coefficient  of  element  n. 

B  Multiple-scattering  amplification  coefficient  of  element  n. 

Zar 
Cnm  Ratio  of  acoustic  impedances,  defined  as:    Cnm  =     "rm  . 


Z 


mm 

ar    .   ^a 


7a!     ,    7a 

Dn  Ratio  of  acoustic  impedances,  defined  as:    Dn  = nnar 

Znn 
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Fa  Scattered  acoustic  pressure  from  a  single  scatterer, 

normalized  to  the  incident  acoustic  field. 

Fb  Scattered  acoustic  pressure  from  a  scatterer, interacting  with 

other  scatterers,  normalized  to  the  incident  acoustic  field 

I  Electrical  current. 

P0  Ambient  pressure  of  environment. 

Sn  Surface  area  of  element  n. 

T  Transduction  coefficient. 

Un  Volume  velocity  of  element  n. 

Un  Volume  velocity  of  element  n,  normalized  to  transducer 

volume  velocity. 

Z  Acoustic  impedance,  pressure  divided  by  volume  velocity. 

Zn  Acoustic  impedance  of  element  n. 

Znn  Acoustic  radiation  impedance  seen  by  element  n  with  all 

other  elements  blocked. 

Znm         Acoustic  transfer  impedance  from  element  m  to  element  n. 

ZfJt  The  total  acoustic  radiation  impedance  seen  by  transducer  in 

the  presence  of  sympathetic  resonators. 

Zao  Transducer  open  circuit  (1=0)  acoustical  impedance. 

^nm  Kronecker  delta  function,  where:    5nm  = 

|  0;  m  *  n  I 

p  Density  of  water. 

H  Ratio  of  the  magnitudes  of  the  scattered  acoustic  pressure  at 

the  surface  of  a  resonant  scatterer  to  the  incident  acoustic 
pressure  in  the  presence  of  other  scatterers. 
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\iq  Ratio  of  the  magnitudes  of  the  scattered  acoustic  pressure  at 

the  surface  of  a  single  resonant  scattererto  the  incident 
acoustic  pressure. 

y  Ratio  of  specific  heats  (air). 

co  Angular  frequency  of  plane  wave. 

coq  Angular  frequency  of  scatterer  resonance. 

C  Defined  as:    £  -  Wf  '  x)  . 

b  ka 

IXnl         A  Nxl    column  matrix,  n  =  1...N. 
[Xnm]       A  NxN  square  matrix,  n,m  =  1...N. 
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